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SUMMARY 


The velocity potential for a triangular wing with subsonic leading 
edges experiencing harmonic deformations in supersonic flow is treated 
herein.- The oscillations considered are such that the amplitude of 
distortion of the wing can he represented by a general quadratic expres- 
sion for a surface. The velocity potential is obtained in the form of 
a power series in terms of the frequency of oscillation. Although only 
terms appropriate for expressing the potential to the third power of the 
frequency are presented, additional terms may be obtained if they are 
desired. The material constitutes an extension of work given in NACA 
Report 1099. 


INTRODUCTION 


Designers of supersonic aircraft are leaning more and more toward 
the use of triangular plan forms for wings and control surfaces. In 
order to obtain information concerning the aeroelastic properties of 
such wings, knowledge is needed of the air forces that may act on them. 

In flutter calculations for wings whose deformations can be calcu- 
lated satisfactorily by simple-beam theory, the use of aerodynamic 
coefficients based on harmonic pitching and translation of representa- 
tive sections of nondistorting, or rigid, wings has proved reasonably 
satisfactory. Simple-beam theory, however, does not as readily apply 
to triangular wings or to swept wings of small aspect ratio and for a 
structural or a flutter analysis of these wings recourse must be had 
to a more appropriate theory. The structural analysis would define 
the structural deformations or the natural modes of the wing, and the 
flutter analysis would involve the use of aerodynamic coefficients 
related to these structural deformations or modes . 
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In references 1 and 2 use is made of classical plate theory to 
develop a method for calculating stresses and deformations of thin, 
homogeneous, cantilever wings of arbitrary plan form. Results of this 
method when applied to a triangular plan form indicate that aerodynamic 
coefficients associated with sinusoidal distortions including deforma- 
tions of the second degree in bending and camber might he useful in 
flutter analyses of triangular wings, particularly if a modal-analysis 
approach is used. 

In references 3 and I 4 - aerodynamic forces and, moments on rigid tri- 
angular wings harmonically oscillating in pitch and in vertical transla- 
tion in supersonic flow are derived by a method of expanding the associ- 
ated velocity potentials in terms of the frequency of oscillation. This 
method of expansion can also he used to derive the potentials associated 
with harmonically distorting wings, provided the form of distortion is 
known. This method, of course, applies to triangular wings with both 
supersonic and subsonic leading edges, hut only the subsonic-leading-edge 
case, which is theoretically the more difficult one, is treated herein. 

The main purpose of the present paper is, therefore, to obtain the 
expanded velocity potential for a triangular wing with subsonic leading 
edges undergoing general second-degree forms of harmonic distortion in 
both spanwise and chordwise directions and to present the first few 
terms of the expansion. Such a potential can be made to yield not only 
flutter coefficients for a distorting triangular wing but also certain 
steady-state and time-dependent stability derivatives. 

For triangular wings with supersonic leading edges the treatment 
is a straightforward procedure, since the boundary- value problem for 
the velocity potential in this case can be satisfied, as shown in' refer- 
ence 5, by simple distributions of sources with local strength propor- 
tional to local downwash. Although this case is not directly dealt with 
herein, the "sonic case," or the case where the Mach lines from the 
vertex of the triangular wing coincide with the wing leading edges, 
which can be considered as belonging to either the case of subsonic or 
supersonic leading edges, is used as a check, in the limiting case, on 
the results presented herein. 

It is of interest to note that the circular plan form in incom- 
pressible flow has been treated for the same forms of harmonic distor- 
tions considered herein (ref. 6). A useful byproduct of such treat- 
ments is that they furnish means by which various strip-theory approaches 
for obtaining aerodynamic forces and moments for a deforming wing, as 
commonly used in flutter analyses, may be evaluated. 
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SYMBOLS 


A,B,C,D,E,F 

A i> B i;> c i> B i> E i.> E i 

A i.> B i.> c i> B i> E i.> E i 


a,b,c,d,e,f 


a i-’ B i-> c i> <p L-’ e i> p i 


t I ! I Jr. 

a i >°± ■ > c l ; e i }*-± 

n , 11 11 , 11 11 _ 

a i > B i j c 1 ; a i ; e i > 



a nm 


2b 


c 


F 


1 


, E 


1 


k 

M 


Ap 

R 


r 

t 


constants defined in equation (4) 

constants depending on Pq defined in' appendix 

coefficients used in equation ( 3 ) to define 
displacement of wing 

constants depending on Pq and M 

constants depending on pg 

functions of 0 , x, and M 
root chord of wing 
velocity of sound 

complete elliptic integrals of first and 
second kinds, respectively, with modulus 

doublet distribution functions 
reduced frequency, boi/v 
free- stream Mach number, V/c 
local pressure difference 
distance, ^(x - - p2(y - 

region of integration (see fig. l) 
time 


V 


free-stream velocity 


W(x,y,t) 


vertical velocity at surface of wing 


W: 


W P 

n,m>"n,m 


constants associated with doublet distribution 
functions depending on Pg and p 
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x,y,z 

Zm(x,y,t) 

p = /# - 1 
A 

I* 1 ! 

P 

P 0 =PA 

0 

CD 

a> = M^ro/vp 2 


rectangular coordinates attacked to wing 
moving in negative x-dlrection 

sinusoidal displacement of a point on wing 


tangent of half -apex angle of triangular 
wing 

rectangular coordinates used to represent 
space location of doublets in xy-plane 

density 

disturbance-velocity potential 
frequency of oscillations 


ANALYSIS 

The Boundary-Value Problem 


The differential equation, referred to a rectangular coordinate 
system moving forward at uniform speed V in the negative x-dlrection, 
that must he satisfied by the velocity potential $ is 






( 1 ) 


The main governing boundary condition, that of tangential flow at the 
wing surface, is 



= W(x,y,t) 



J^Zm(x,y,t) 


( 2 ) 
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where W(x,y,t) and Zm(x,y,t) represent, respectively, the vertical 
velocity (or dovnvash) and vertical displacement of any point (x,y,z) 
of the wing. Let the amplitude of vertical displacement of the wing he 
determined hy a general second-degree equation in . x and y; then the 
displacement of any point of the harmonically oscillating wing may he 
expressed as 


Zm = e^^ax 2 + hy 2 + cxy + dx + ey + f ) ( 3 ) 

whfere the coefficients a, h, c, . . . may he considered as complex 
quantities in order to permit the inclusion of phase differences between 
the different components of the displacement . By substituting the 
expression for Z m of equation ( 3 ) into equation (2), an expression 
for the dovnvash is obtained which may he written as 



W(x,y,t) 


= e^^Ax 2 + Bp^y 2 + Cpxy + Bx + E[3y + F) (h) 


The commonly used Mach number relation pi s Jlfe- - 1 has been introduced 
into the coefficients of the y-terms in equation (If) for convenience in 
the subsequent analysis. The coefficients A, B, C, . . . are related 
to the coefficients a, b, c, . . . in equation (3) in the following 
manner; 


A = ina 



E = 


Vc + itoe 

P 



D = 2aV + ind F = Yd + im£ 


Of course, for some purposes it may be more logical to start with a 
prescribed dovnvash or vertical velocity rather t h a n a prescribed 
deformation. 

As is usual in dealing with linearized aerodynamic problems, the 
velocity is conveniently expressed as a sum of separate effects asso- 
ciated with the different terms of the dovnvash expression, ' equation (If), 
namely 
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0“0A + 0B + 0C + $D + $E + $F 


( 5 ) 


With the potential expressed in this form attention may he directed to 
the derivation of these subsidiary potentials $A; $Bj> and so forth, 
individually . 


Derivation of Velocity Potential $ 

Method of solution .- The boundary -value problem for each of the 
subsidiary potentials in equation ( 5 ) can be satisfied by the method 
discussed in reference 4 of expanding the velocity potential to any 

given power v of the frequency parameter d> = • The procedure 

for obtaining any one of these potentials is essentially the same aB 
that for obtaining any other. It will therefore suffice to discuss the 
derivation of only one subsidiaiy potential, say 0A* ^h e observations, 
pertaining to pulsating doublets, made in the following paragraph may be 
helpful in the subsequent outline and discussion of the derivation. 

When the potential of a doublet that is pulsating with frequency tn 
and moving at a uniform speed V is expanded in terms of 00 , certain 
terms of the expansion are found to contain as a factor the so-called 
steady-state doublet potential z/R^ (hereinafter referred to as sin- 
gular terms ) , whereas all other terms contain factors of the type 

(a/M) 2m- ^zR^ m- 5, m ^ 2 (hereinafter referred to as nonsingular terms). 
Of these terms, only the singular ones contribute directly to the poten- 
tial for the airfoil; nevertheless, both singular and nonsingular terms 
are, in general, required to satisfy the differential equation, equa- 
tion (l), to a given power of 55. Furthermore , when these expanded 
potentials are employed in satisfying the boundary conditions, the non- 
singular terms give rise to downwash that, as indicated by the fac- 
tor (cs/m)^ 111- ^, is of higher degree in d> than the downwash arising 
from the singular terms, an observation that provided a key to the step- 
by-step procedure adopted in reference 

In this procedure the singular terms can be grouped and weighted 
so that they alone satisfy the specified downwash conditions. At the 
same time the downwash arising from the nonsingular terms is to be 
canceled to the required power of d) by the introduction of other 
doublet potentials in order to preserve the downwash conditions sat- 
isfied by the singular terms . It is through this canceling process 
that the nons ingular terms give rise, indirectly, to any contributions 
to the potential for the airfoil. 
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The distribution functions .- In accordance with the preceding dis- 
cussion, the potential under discussion can conveniently he considered 
as a sum, namely 


" $1 + 01 + $1 + • • • 


( 6 ) 


where $1 60 treated that the singular terms contained therein_alone 

can be made to satisfy the boundary conditions and the terms Sj_, . . 

are successively required to cancel the downwash arising from the non- 
singular terms of $1 • The term , which may be considered as a 

summation of weighted doublets expanded to the vth power of o> nnd 
integrated over the appropriate portion r (see fig. l) of the wing 
plan form, may be written as 



= Ae imt a nl x n G n (x,y) (7) 

n=0 


where 

R = |/(x - i) 2 - p 2 (y - t 1) 2 - p2z 2 

and 


s-nl 



(-idk) v - n 
+ (v - n)’ 


( 8 ) 


(-1) 131 - 1 

(2m - 2)1 




icux + . 


(-tox)^ 2 - 11 - 2111 
( v + 2 - n - 2m) ! 


a nm — 


(9) 
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The functions G n (£,r]) represent weighting or distribution functions 
that are to he determined so that the first set of integrals in the 
brackets in equation (7) (which, as may be noted, employ as kernels the 
singular terms of the expanded doublet potential) satisfy the prescribed 
downwash condition (see ref. 4 for a more detailed discussion), namely 



( 10 ) 


Corresponding expressions for (f> $j_, . . . may be written as follows: 


~7 Ae 


iorfc 


lim 


8/S 2 


v-2 


rt z-e>0 MSI 1 ! 2 / n^ 


Kil ff ®n( ^ d n + 

r 


2v-l+(-l) 1 


m=2. 


G n (i > 'n)| n+3 R 2m " 3 d| drj 


n=0 2Mr 


(H) 




icut 


lim 


ap-\ 

0 Hwy 


n=0 


a nl JJ G n (|,T|)£^-? d l dl l + 


2v-3+(-l) 


m=2 


> a nm / / G n (^T 1 ) l n+ 5 R 2m-3 d£ dtj 


= Ae 


laA 51 a nl xn+i ^Gn( x ^y ) 

n^ 


( 12 ) 
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■where the functions Gn(|,i]), G n (£,n), • • • are to he determined so 
that 0 L , . . . yield a downwash that -will cancel the downwash 

arising from the second set of integrals in the equation for 0-j_, equa- 
tion (7), which employ nonsingular terms as keraei&. ■ 

Since the following indentity holds (see eq_. (13) of ref. 4 ), 

. v 

21 xIla nl = 1 
n=0 

it follows that the integral equations which determine successively 
the distribution functions Gn» Gn, Gn, • • . (n = 0, 1, 2, 3, . . .) 
so that' the condition of tangential flow for is satisfied are as 

follows : 



ff Gn(i>ri)£— ^ d| dri = Ax 2 (x n ) = Ax n+2 
r 


(13) 



► n+3 

dg dr, = 


.lim ~ [f G n (|,T])| n+1 R dg dT) 
Z — ->0 bz*# 


(14) 


^2 t n +5 ^2 

Hm — - // G n (|,ri)l— d| drj = lim — 
z — K) dz^ r R z — >0 dz c 


iff Gn(l,n)£* +1 R 3 <^ dTi 


K 


ff G n (£; T l)£ n+5 R <11 IT) 


(15) 


These integral equations may be solved by the method discussed in the 
appendix of reference 4 where, in view of integration difficulties, it 
is found necessary to consider not only the downwash that is to be 
satisfied but also the derivatives of the downwash. This method, 
though simple in principle, is rather lengthy to apply; consequently, 
only the results of the solutions will be given herein. Furthermore, 
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the potentials •will hereinafter he restricted to terms up to the third 
power of c5_ (that is, v = 3)- To this power of B the distribution 
functions G n , equation ( 15 ), are not required. 


Expressions for Gn and Gn are found to he as follows; 


G 0 = I 2 


aq + a x JliUA 2 ! 2 - -n 2 


Gi = i 2 


G2 = I 2 


g 3 = r 



^7 + A8 + *9 ^^)A 2 i 2 - t ) 2 


> (16) 


% = & 2 


% * A! + As E^V 2 6 2 - I 2 
v r r / 


Gl = s 2 


% * b - 1,2 ^ 


+ a 5 - 1‘ 


The coefficients Ap and Ap in these equations are functions of 
Pq = (3A only, where A is the tangent of the half-apex angle of the 
triangle (see fig. l) . These coefficients arise in the process of 
solving the integral equations and are defined by systems of linear 
algebraic equations which are given in the appendix. • 


The integral equations for the distribution functions associated 
with 0g, 0 q, • • • , equation ( 5 )* are of the same general form as 

those associated with 0^. The main difference is that the term on 
the right side of equation (13) is replaced by the term appropriate to 
the downwash under consideration. The solutions yield the following 
distribution functions: 
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For 


for 0 C : 


GO = I 2 


G X = | 2 


G 2 = | 2 


G3 = | 2 


G o “ £ 2 


Gi = £ 2 


+ Bi ££W 2 E 2 - ^ 


Bq + B X 


+ B 3 ^|^)/l 2 i 2 - T1 2 

+ ®5 ^2~ + B 6 ^)A 2 £ 2 - Tl 2 
^7 + B8 ^2“ + B^ ^^^ A 2 ^ 2 ~ T l 2 
^0 + Bi 2?a£ + Bs ^)A 2 ^ 2 - n 2 


A 


B 3 + b 4 ^ + B 3 ^)A 2 fi 2 - n 2 

^ 1 r 


Gq = m (c 0 n|a 2 ^ 2 - n 2 ) 


G 1 - Ptj^ 


g 2 - 


g 3 = PTJl 


J 01 ♦ Cg 


(03 + ci(. ^)^T 7 

(«, + c 6 ^ + c 7 


(17) 


> (18) 


Gq ~ P 1 !! 


Co + Cl - n 2 
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and, finally, for 9 ^.: 

g 0 « (f 0 4 2 ! 2 - n 2 ) 



g 2 = 

g 3 = 

Gq = 
Gl = 



+ f 3 


fa + *5 




+ Fj 



(21) 


The coefficients BjL> Bj_, Cj_, Cj_, . . . are functions of Pq = pA. 
They are defined by linear simultaneous equations given in the appendix 
following those defining A^, A^. 


Expressions for the potential .- With the use of the distribution 
functions, the potentials can be written by inspection and .comparison. 
For example, from equations ( 7 ) and (ll) the expression for $a in- 
terms of Gn and Gq is 


0A = Ae ia>t ]>2 a^x^^y) + Ae^ 
n=0 


2Z ~~n a nlx n+ 2 G n (x,y ) 
n=0 2Mr 


( 22 ) 
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The potentials $B> $C> • • • are represented by similar expressions. 

After terms have been combined with respect to tu, the subsidiary poten- 
tials to the third power of to can be written in the following forms: 

$A = Ae” 5 ^ ^/a 2 x 2 - y2 |^x 2 + a^p 2 y 2 - im {a^x? + a^xp^y 2 ) - 

B 2 (a.]^ + a 5 x 2 g 2 y 2 + agpV - ) + iu?(a^x5 + agx3p 2 y 2 + a^xpV^) (23) 

= Be lcut \/a 2 x 2 - y 2 |bpx 2 + Bpp 2 ^ - ico (b 2 x3 + bjxP 2 y' 2 y - 
^(bjj.x 4 ' + b 5 x 2 p2y 2 + bgpV*) + io>3(b^x5 + bgx^pSy 2 + b^xpV 4 ") (24) 

$C = Cye lxut ^A^ 2 - y 2 jcgx - iiu^c^ 2 + c 2 p 2 y 2 ^) - to 2 ^^ + c^xp^ 2 ) + 
iio^ (c^x^ + cgx 2 p 2 y 2 + c^pV* - ) (25) 

= Be iait \/a 2 x 2 - y 2 ^x - ico^d^x 2 + d 2 P 2 y 2 ) - to 2 (cU,x3 + dj ( xp 2 y 2 ) + 

^ ( d 5 x ^ + <%x 2 p^y 2 + (26) 

(27) 

(28) 


$E = Eye 1 ^ /a 2 * 2 - y 2 jep - itoe-jx - to 2 ^^ 2 + ejP 2 y 2 ^ + 
ioj3^x3 + e^xp 2 y 2 yj 

$F = Fe icot \/a 2 x 2 - y 2 Jf 0 - itofpx - a5 2 (f 2 x 2 + f^P 2 y 2 ) + 
ia>3 (^f 4x3 + f^xp 2 y 2 yj 
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A-i; 


In these equations ap, bp, . . . are linear combinations of 
Ap, Bp, Bp, • . . and can he further expressed as 


a i = a i' 


M 2 


hi = bp’ 



Expressions for the quantities ap’, ap" , bp', bp", . . . are given in 
the appendix in terms of the quantities Ai, Ai, Bp, Bp, . . . . In 
addition, values for these quantities for particular values of pQ 
throughout the range 0 ^ Pq % 1 are given in table I, and, for con- 
venience in interpolation, they are shown plotted in figure 2. Thus, 
for a given triangular wing and a given Mach number, such that 

Pq 2 = p 2 A 2 < 1, the coefficients ap, bp, . . . fp in equations (23) 

to (28) can be assigned their numerical values with the use of table I 
and figure 2. 

It may be noted from the boundary condition, equation (A) , that 
the subsidiary potentials fa, fa, and fa are associated with motions 

of rigid wings; and that the potentials fa and fa, equations (26) 
and (28), could be obtained directly from reference ^ by comparison 
with the potentials denoted therein by fa and 0 a , respectively. Also 
the parameter eg in equation (27) agrees, as it should, with the param- 
eter associated with constant rolling motion of a triangular wing (see, 
for example, ref. j). 

The sonic case .- At f3A = 1 or A = l/p, which is the condition 
at which the Mach lines from the apex of the triangle coincide with the 
leading edges of the triangle, equations (23) to (28) reduce, respectively, 
to 
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Wb) 3A =i 


^c) pA =l 


2A\6c2 - p%2 oiojfc 
P* 




^-(llx 2 + p2y2) - (l7x3 - 5xp 2 y 2 ) - 


S 2 


1575M 2 


io^ 


155925M £ 


|0-7 + 29M 2 )x^ - (ik + 13M 2 )x 2 p % 2 - (3 - 4#)pV^ 


[(759 + 4oiM 2 )x5 - (858 + 137M 2 )x3p2 y 2 


(99 + 36M 2 )xpvJ I 


( 29 ) 


^ •** {i §(* 2 * ^% 2 ) - IfC - 5 - »*0 - 


m 2 


1575M 2 


ia^ 


31185 M 2 


|(3 + 7M 2 )x^ + (l k + M 2 )x 2 p 2 y 2 _ (17 _ 12M 2 )pV] 
[(33 + 23M 2 ) X 5 + (66 - 23M 2 ) x 3p 2 y 2 _ 


(99 - 60M2)xpV 1 ' 


(50) 


2Cyy£2 - p2 y 2 J 

P* 


^ . p 2 y 2) 

15 315 V ' 


315M 2 


f?l + M 2 )x3 - xpV] + . . ^ . 1 ^ — 1?66 + l4M 2 )x^ - 
Li J 51975 M 2 L 


— V 

(77 - i7M 2 )x 2 p 2 y 2 + (11 - 6m 2 )pV^_ 


(31) 
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WpA=l 


2D 


\)x 2 - p 2 y 2 




P* 


# - frK 4x2 - P 2 ^ 2 ) - ———2 |C T + HM 2 )^ - 

3 45 XT rm 2 [A 


315M C 


(7 + 2M 2 )xp 2 y 2 ~] + — ^-|?54 + 26 m 2 )x^ - 
4725M 2 L3 


(63 - 3M 2 )x 2 p 2 y 2 + (9 - 4M 2 )pV- 


(32) 


(^) = e i*t . 2^ _ S |( 5M 2) X 2 _ 

pA" 1 prt 1 3 15 rtrm 2 ,v ' 


315M'= 


(j - 4t^)p 2 y 2_ ] + ■ (9 + M 2 )x3 - (q - kM^xp^ 2 

- 1 9A5M 2 L 


(33) 


(9(f) 


= e tot f x . . js1|7 5 + ^2 . 

_1 P« 3 90M 2 >- 

(5 - 2M 2 ) P% 2 J + - la ^ g (7 + 3M 2 )x3 - (7 - 2M 2 )xp 2 y 2 jj* (34) 


Equations (29) to (34) have also "been obtained by integrating the 
expanded potentials of appropriate sources over the region occupied by 
the wing for this case of pA = 1. They therefore serve as a check on 
the results given in equations (23) to (28). 

Introduction of the reduced- frequency parameter into the expressions 
for the potentials . - In applications it is usually desirable to refer all 
lengths to some convenient reference length on the wing and to introduce 
a reduced-frequency parameter. For this reason it may be desirable to 
write the potentials in a slightly modified form. If, in the present 
case, the root chord 2b is chosen as a reference length and the varia- 
bles x -and y in the potentials are employed in a new sense to mean 
that they are referred to this length, the expressions for the potentials 
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can be written in terms of the reduced- frequency parameter k = bcs/v 
as follows: 


fa = 2b(Vb 2 A)e :ix0 ^ \/a 2 x 2 - y 2 


a^x 2 + a 1 p 2 y 2 - 21 - M . k ( a 2 x^ + a^xp^ 2 ) 


3‘ 


+ a 5 x 2 p 2 y 2 + agpVO + 


8 ~ ^ k + a 8 x5 P 2 y 2 + ^pVO 


(35) 


0 B = 2b(kb 2 B)e ia3t i/a^ 2 - y 2 
kM^k 2 


b 0 x 2 + b B p 2 y 2 - 21M k (bpx3 + bjxp 2 y 2 ) 


-(b^x^ + b-x 2 p 2 y 2 + bgpV^) + 


8iM 6 k3 


-^—(brjx5 + b Q x3p2y 2 + b^xpVO 


(36) 


fa = 2b(kb 2 C)ye lxut /a 2 x 2 - y 2 


■V - ^(?1* 2 + =2l3 2 y 2 ) 


—^(05x5 ci^y-P^-y 2 ) + (c^x^ + cgx 2 p% 2 + CyP^A) 


(37) 


fa = 2b(2bD)e layfc /a^ 2 - y 2 


*0* - ^^( d l x2 + d 2P% 2 ) " 


^ k - - (d^x3 + d^xp 2 y 2 ) + 8lM -^ (a^x 1 *- + dgx 2 p 2 y 2 + d^pV - ) 


(38) 
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(fhg = SbCabEjye^ \/a 2 x 2 - y 2 


2iM 2 k 4M^k 2 

•" ' IT" eiX - — 


(egX 2 + e 5 p%2^) 


8iM^k^ 


(ejj.x3 + e^xp% 2 ) 


(39) 


0P = 2b(F)e tot A 2 x 2 - y 2 


„ 2iM 2 k . W*-k 2 / P _ ? \ 

f 0 — fyX — (f^ 2 + f 3P y J + 


8iM^k3 


(f^x3 + f 5 x^ 2 y 2 ) 


(40) 


(it is to tie noted in these equations that the expressions within the 
first parentheses in each equation, namely those containing the coeffi- 
cients A, B, C, . • . , are expressions for downwash and hence possess 
the dimensions of velocity.) 


REMARKS AND DISCUSSION 


The local force (positive downward) or pressure difference Zp 
between the upper and lower surface at any point (x,y) on the wing is, 
in terms of nondimensional coordinates, 



(4-la) 


or, for the oscillating case under consideration, 


Ap = - 


pv/a 

b \^x 


2ikW 


(4-lb) 
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By substituting tbe potentials given by equations ( 35 ) to (40) into 
these expressions , the pressure distributions can then be found for the 
types of motions or downwash considered in this paper. It should be 
pointed out that the expressions for potentials developed herein for a 
triangular wing can be Used to calculate local pressures oh certain 
other plan forms, as indicated in reference k, which can be formed from 
the tri angular wing by cutting the trailing edges in such a way that 
they lie ahead of Mach lines emanating from all points of the trailing 
edges . 

The pressure distributions obtained in this way can then be inte- 
grated to obtain the desired spanwise distributions of the lift and 
moment, as well as any integrated quantities which may be of interest 
in the structural analysis of the given wing. For instance, in a modal 
analysis, for which the results of this paper would be particularly 
useful, the aerodynamic forces usually enter in the form of integrals 
over the wing of the product of local pressure difference and the local 
distortion in a given mode, these integrals representing virtual work 
done by the aerod ynami c forces. The virtual-work integrals arising 
from ahy or al 1 of the potentials presented in the foregoing analysis 
can be easily reduced to involve only integrations of functions of the 

form x n (x2 - a^)* 1 / 2 and can therefore be readily evaluated. It is 
recognized that, if the wing root is considered to be rigidly fixed, 
some of the displacement terms considered herein will not satisfy 
boundary conditions at the wing root for either plate or beam theory. 

If, however, such terms are useful in approximating a known mode, shape, 
these types of disparities may be overlooked with the interpretation 
that they imply large stresses at the wing root especially since with 
triangular wings such stresses are known to exist. 

Inasmuch as the types of integrals required vary for different 
methods of aeroelastic analysis (depending on whether the theory for a' 
simple beam or that for a more refined representative structure is used) 
and, inasmuch as the limits of these integrals depend on the exact plan 
form, further manip ula tions employing equations .(^-l) to obtain some form 
of force and moment coefficients are not presented herein. However, in 
order to illustrate the results obtainable by the method of this paper, 
the spanwise variations of the unsteady-lift derivatives associated with 
parabolic bending (obtained from $ 5 ) ara shown in figure 3 for a 
lt-5° delta wing undergoing parabolic bending at k = 0.1 and at a Mach 
number of 1.2. These results may be used to give some indication of 
the way coefficients based on the true downwash conditions for the 
distorting wing compare with those obtainable by a simple strip-theory 
analysis based on the coefficients associated with rigid-wing motions. 

In order to effect such a comparison, the results of multiplying the 
lift derivatives associated with translation (obtained from $F ) by the 
ordinates of the assumed parabolic mode shape (y^) are also shown in 



NACA TN 3009 


21 


figure 3 for the same value of k and M. The main features to note 
in this comparison are the differences in magnitude of lift at the -wing 
root and at the points of maximum lift. These differences are, of- 
course, associated with induced effects due to bending. In order to 
determine the nature of these induced effects on calculated flutter 
speeds, further investigations are required. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June k-, 1953 * 
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APPENDIX 

COEFFICIENTS ASSOCIATED WITH THE VELOCITY POTENTIAL 


This appendix contains definitions of the coefficients that have 
arisen in the foregoing analysis . 

The coefficients denoted by A±, Ap, B]_, Bj_, . . . in equations (16) 
to (2l) of the body of this paper are defined by systems of linear alge- 
braic equations which arise in the process employed for solving the 
integral equations involved. These algebraic equations are listed sub- 
sequently; however, they in turn involve- a set of cumbersome coefficients 

that are denoted by W^ m and W^ m . These coefficients are defined as 
follows : 

Application of the transformations (see appendix of ref. 4) 


(3x(e - a) 

N = 


1 - ( 3 ^ 0 ^ 


q. = 


x(l - p 2qe) 
1 - { 3 ^ 0 ^ 


0 


_ y 


X 


0 


n 


i 


to the integral equations (and certain derivatives thereof) which define 
the distribution functions (see eqs . (13) and (lh) for $a) results in 
the definitions 
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W 


P 

n,m 


(n + 2)p 2 r +A 0 m \/A2 - do 

rtx 11 J-K .(1 _ g2g-2)3/2 


<? 

11m - — 

0— >0 de p 




Nn+l 
cosh t ) 


cosh r dr 


W P — P 

w n,m 


X cM 2 - o 2 


dcr lim 




itx n+1 /l _ r 2 ct 2 9—^0 S0 P ^0 


^ cosh' 


-11 


IT , >T , . n+l r 

( q - N cosh t) 2q 


(n + 4)N cosh -rjslnh^T dT 


(Unless m and p are either both even or both odd, w|^ m and $P m 
vanish owing to symmetry . ) 

With F 1 and E 1 denoting complete elliptic in tegrals of the 
first and second kinds, respectively, with modulus '/l - p Q 2 , explicit 
values of W p ^ m and wP^ m required to write the equations for Aq, Aq, 
Bf, Bf, . . . are 



= E' 







2k 
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K 1 - p o 2 ) : 


2 11 


(5P 0 2 - 5p 0 4 >’ 


+ 


(2 - 10P 0 2 




»il - 1 jftpo 2 + Po’ t >' + ( 2 - 2P 0 2 + 2Po 4 ) eT 1 

' (i-Po 2 ) L J 







M°,0 * RtPq 2 - SXPq^ + 12P 0 6 >' + (6 - 55P 0 2 + 

6(1 - Po 2) 3 L 

W - * p 0 *>] 


— — ? ^P 0 2 - Pq^f' - (5 + 7p 0 2 - 2P 0 1| ') E ] 

C 1 - p o 2 ) 


w| A - 1 Jpo 6 - 9P 0 k - 2Pq 2 )f 1 + (k - 5Pq2 + 15Po k - SPo^j 

2(l - P 0 2 ) 
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^[(po 2 - 9P 0 k y - (2 - 7P 0 2 - 3P 0 4 > 

- Po 2 ) 







+ 9P 0 4 


- 3p 0 2 )f' + (6 - I5p 0 2 + 5p 0 I <- - kp 0 6 )E ' 







3lP Q k + 21 p q 6 )f' + (2k - 65 p q 2 + 



<0 


8 ( 1 - eo 2 ) 


(56p 0 2 - 92p 0 j + + 72 p 0 6 - 20p 0 8)f’ + (8 - H7P 0 2 + 


202p 0 ^ - 


1^9P 0 6 


itOpQ 8 


5 )e' 



— - — rlO-TPo 2 - V + Po 6 >' - ( 4 + 18 P 0 2 - V + 2 p 0 6 >' 
<1 - Po 2 ) J 
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w 4,o — r[( 8p o 2 + 8p o 4 ) F ' “ ( x + + p o 4 >’ 

(-Po 2 )^ J 

wj 1 = -(2p 0 2 + 24 P( ^ - l4p 0 6 + 4 p 0 8 )f ’ + (4 - 5p 0 2 + 

2 ( 1 - Po 2 ) 

38p 0 ^ - 29 p 0 6 + 8p 0 8 )e ' 



W4 2 - T- -(2P 0 2 - 10p 0 h - 10p 0 6 + 2 p q 8 )f ' + (4 - 13P 0 2 + 

2 (1 - p 0 2 /L 

2p 0 h - 13P 0 6 + 4p 0 S )e 1 



2p 0 ^ + 17 p 0 6 ) f ’ 


- (2 - 8p 0 2 + 
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Wj 


4,3 


-(l2p 0 2 - 42p Q 4 + 72 Pq 6 + 6 Pq 8 )f ' + (24 - 87 p 0 2 + 

C 1 - O L 

ll4p 0 4 - 15P 0 6 + 12p 0 8 )e’J 


TT° - P 0 t A 

M ““8 

240 


W 


4,4 


23 2 (l - P 0 2 ) 
4 p 0 8 )e] 


( pQ k - 2P 0 6 + 17P 0 8 )f' - (2p 0 2 - 8p 0 ^ + l8p 0 6 + 


4,4 - 5 


0 - Po 2 ) 


(56p 0 8 - 92P 0 6 + 72p 0 4 - 20p 0 2 )f ' + (40 - !49p 0 2 + 


202p 0 ^" - U7P 0 6 + 8p 0 8 )e ’ 


yO — 1 


4o(l - p 0 2> 


(4O0p 0 2 - 833 p 0 4 + 994p 0 6 - 553 p q 8 + 120p 0 10 )e’ + 


w , 


5,0 


(4o - 859 p 0 2 + 1910 Pq^ - 2115 p 0 8 + U36 p 0 8 - 240p 0 10 )e’ 


(li5P 0 2 - 2p 0 ^ + 19 p 0 6 - 4p 0 8 )f 1 - (20 + 151p 0 2 - 


2 (l - Po 2 ) 5 L 
74p 0 ^ + 39p 0 6 - 8 p 0 8 )e J 
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«5,0 " 


M 


h n 


C 1 -^ 2 ) 5 


(55p 0 2 + 7^p 0 4 - Pq 6 )?' - (5 + 108p q 2 + i T p 0 4 _ 2 Pq 6 )e 1 


Wi 


1 


5,1 / n5 

8(1 - P 0 2 if 


(20 p 0 10 - 9lp 0 8 + 15^P 0 6 - 203p 0 4 - 8 p 0 2 )f ' + 


(16 - l6p 0 2 + 323p 0 4 - 342p 0 6 + i 87 p 0 8 - 1k)p 0 10 )e' 


Wi 


3 _ 3P‘ 

5,1 “ 


?(1 - P 0 2) 


5H 


(i^p 0 2 - 131P 0 4 + 2p 0 6 - 3p 0 8 )f ' - (8 - 53p 0 2 - 


108p 0 ^ + 31P 0 6 - 6 p 0 8 )e" 


Wi 


5,1 


3F 


C 1 - p o 2 )' 


l^o 2 - ^ kp Q k “ 55P($F’ - (2 - 17P 0 2 - W8p 0 k - 5p 0 6 )e] 


w® = Q w 2 

5' 2 aop* 5,° 


Wi 


5,2 


(8p 0 10 - 35P 0 8 + H2P 0 6 + k 9P 0 h - 6 p 0 2 )f' + 


2(1 - P 0 2) 5 

(l2 - 1^7P 0 2 - 22p Q ^ - 127 p q 6 + 72P Q 8 - 16 p 0 10 )e 


w; 


5,2 


3P C 


0 - - po 2 y 


|p 0 2 - 2P 0 ^ + 131P 0 6 - ^Po 8 ^' - ( 6 - 51P 0 2 + 
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w 5,3 


(6p 0 l° - i^9p 0 8 - H2p 0 6 + 35p 0 ^ - 8p q 2 )f ’ + 


K 1 " PQ 2 ) 5 

(16 - 72p 0 2 + 127P 0 ^ + 22p 0 6 + VfP 0 8 - 12p 0 10 )e^] 



w 2 ^ = _ I (3P 0 ^ - 2P 0 8 + I3ip 0 8 - 4 p 0 ^)f' - (6p 0 2 - 

2 P 2 (l - p 0 2 ) 5 

3ip 0 4 + U 08 p 0 6 + 53p 0 8 - 8p 0 10 )e’ 


W5 k = - K8p 0 10 + ^Pn 8 - i^P 0 6 + 9iP 0 ^ - 20 p 0 2 )f' + 

0 - Po 2 ) 5 l ^ 

(%0 - i87p 0 2 + 3 k 2 p 0 k - 323 p 0 6 + i6p 0 8 - i6p 0 10 )e ' 
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l 

6(1 - p 0 2) 2 
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1 
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— — - (7Po 2 + Po^) F ' - ( 2 + 7P 0 2 - Po^E' 


w 


2,0 


- Po 2 > 


w*«- = £ 

"2,0 


C 1 - Po 2 ) 


3.11 


(5P 0 2 + 3p 0 ^>' - (l + 7 p 0 2 >' 


4,1 = P 2 ^ 


3,2 


tj3 _ 

W 2,l 


,2 r- 


(1 - Po 2 ) J ^ 


(p 0 2 - 9P 0 k ) f ' - (2 - YPq 2 - 3pq 4 V 


“2,2 “ 3 [(3Po^ + 5P 0 6 >’ - ( 7 p 0 4 + P 0 6 S') 

24|3 2 (l - P 0 2 ) - 1 


“ 1,2 


K 1 - PO 2 )' 
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C 1 - p 0 2 )' 


“3,0 rl^Po 2 ' ‘‘H 4 + Wp o 6 - 12 p 0 8 ) r ’ - (“tpq 
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126p 0 4 + 91p 0 6 - 24p 0 8 )e ’ 


*3,0 - - ~r^ — -x[v 


6 ( 1 - Po 2 ) 


(33P0 2 + l^Po 4 + Pq 6 )e - ( 6 + ^Pq 2 " ^ 0 ^ + 2 Po 6 )e ' I 
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* 3,0 ‘ <0 


3,1 


r (lTPo 4 “ 2 Po 6 + p 0 8 ) f ' - (^P 0 2 + i8p 0 4 - 8p 0 6 + 2 p 0 8 )e ^1 

8 ( 1 - P 0 2 ) J 


tj 3 , _ P £ 


¥ 


3,1 


tj 5 

"3,1 
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w 3,2 


2(1 - P 0 2 ) 

p^ 
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tj 2 _ 

w 3,2 
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- Mj-Pq 1 *' + - ^Pq®)!'’ - (24 - 




The necessary equations for k±, A.±, B^, IS^, . . . are, in terms of 

™ and- as follows: 

n,m n,m’ 

For $a : 


V»2,0 + A l p2H 2,2 
Aq W 2,0 + A l3 2tf i,2 



A 2 W 3,0 + A 3 p2w 3,2 
A 2 W 3,0 + A 3 p2ff 3,2 



Al^W^o + ^3^2,2 + A 63 
V^O + V^,2 + 

\ w J,o + A 5 p2w i[,2 + a 6 p Ma 



A 7^,0 + A 8 p2w 5,2 + a 9^ 1< ^5,^ 
*7*5,0 + k Q^%2 + k 9^,h 
a 7 ^ 5,0 + k 8 fi 2 ^>,2 + a 93^,4 
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Aq^q + A l ^,2 + A 2 P H,4 = + A 1 ^ 2,2 
V^O + A 1 P% 4,2 + A 2^H,k = Vl,0 + A lP^2,2 
AO^,0 + A l p2w i,2 + A 2PH,4 = Ao^ 2,0 + A l p2 ^2,2 

a 3^,0 + A 4 p2H 5,2 + A ^^5,b = A 2^3,0 + A 3 p2 ^3,2 
Aflo + A 4 P%| j2 + A^^k = A 2*3,0 + A 3^,2 
^ 3 ^ 5,0 + A 4 p2 ^,2 + ^ p 4 w 5,4 = a 2^3,0 + A 3 p 2 ^5,2_ 
for <j 6 B : 


> 


> 


V*2,0 + B l p2w 2,2 “ 0 

B 0 W I,0 + B 1 P% 1,2 * 2p2 . 

B 2^3,0 + B 3 p% 5,2 = 0 


b 2 W 3,0 + B 3 p% 3,2 - 2p2 


B^wj^o + 4,2 + B 6 p w 4,4 “ 0 

Bi^o + B 5 p2w 4,2 + B 6 p MU = 2p ‘ 

b 4 w 4,0 + B 3 p ^,2 + B 6 p H,4 = 0 


> 
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b 7 w 5,0 + 



* 7 * 5,0 + 

8^,2 

+ BgiA^,, 

+ 

O 

^tlT\ 

lb 

b 8 b X>2 

+ 

r^U 

B 0 <0 + 

B 1 B X,2 

+ B 2 

B X,4 

+ 


- 4 2 
+ b 2 p w 4A 

vto + 

VX,2 

+ B 2 

sH,k 

%«5,o + 

V%° )2 

+ % 

B X,4 

%«5,o + 

%P 2 w§,2 

+ Bcj 

B X,4 

% w 5,o + 

¥^,2 

+ B^ 

b X,4 

for (f > c : 




C 0 W 2, 1 = 

1 



<>1*3,1 + 

c 2 p 2 w 3,3 

= 1 





► 

°l w 3,l + 

02» 2 ^,3 

= 0 

J 

1 

C 3 W 4,1 + 

VH,3 

= 1 

l 

c 3^,l + 

c 4 b * 4,3 

= 0 

-> 

r 


2p‘ 


> 


0 




B c$U + 

¥2,0 + B lP^2,2 J> 

B 0^2,0 + B l p2 ^2,2 

B 2^,0 + B 3 P ^ 3,2 
B 2^3,0 + B 3 P ^3,2 
B 2»3,0 + * 3^,2 


> 
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°5 W 5A + C 6» % 5,3 + ‘VMa 
°5 W 5,1 + C 6»S,i + vHs 

°5 W 5,1 + c 6P 2 ' , 5,3 + c 7^ w 5,5 



* c lf>%3 - c 0 w 2,l 

S w 4,l + = c cp2,l 




C 2 W 5A + ¥ % 5,3 + 

®3 p ^»3 + ®if pltw 5)5 


C 2 W 5,1 + 


C 1«3,I + °2P^,J 
c l^,x + ^P^A 

^2 W 5,1 + °3P 2 ' , 5,3 * ^^ 5,5 = c l«3,l + C 2P^3,3 
for (ft-Qi 


> 


D 0 tJ I,0 


1 


D 1 W 2,0 + D 2P^ i 2,2 

D 1 W 1,0 + ^ 2 «2,2 



d 3 w° ) 0 + V^.a 

d 3 w 3,o + "A,2 


1 


> 


0 
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D 5^,0 + * 6^,2 + VMU 

d 5 w ?,0 + d 6P 2 w 5 , 2 + 
d 5<0 + D 6P^it,2 + ^P^^ 

V^o + = V?,o 

> 

^,0 + h ^,2 = V ?,0 



»2<,0 + ^^,2 + ^ pM ^ I *. 

*A,0 + ® 3p 2w f,2 + ^- P ^^ it - 

® 2<0 + + vNu 


D 1 ^ 2,0 + D 2 p2 ^ 2,2 
D 1 ^ 2,0 + *2?% 2 
D l * 2,0 + D 2P ^ 2,2 


> 


for 0 ^: 


E o w i,i 


E 1 W 2,1 


E 2 * 3,l 


E 2 « 3,l 


= 1 
= 1 

+ ^3,3 


:: 


E lX,l + E 5 p2 ” 4,3 “ 1 

V'il + e 5 P % 4 ;3 = 0 
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The coefficients a i' * a.^" , bj 1 , b.^", . . . used in the expres 

sions for the velocityjpotential (eqs. (23) to (28)) are defined in 
terms of A-j^ A^_, B-j_, B-^, . . . as follows: 

a o' _ Ao 

II _ 

ao = 0 

a l ' = A 1 

II _ 

a^_ =0 

&2 ' = Ao - A2 

II ^ 

&2 =0 

a 3 ' = A 1 “ a 3 

II _ 

=0 

. ^ A . 

& k _ _ a 2 + _ 

V = - | Ao 

1 A 1 A 3 

a 5 - T - A 3 + 2 

a 5 " = - | A-l 

. a 6 
a 6 - T 

11 1 r 

a 6 = " g A 2 

, Ag A 2 A]j a y 

a 7 "T“T + T“T 

a Y = - |(aq - A3) 

a 8 ' = _ ^5 + h _ ^ 

6 2 2 6 

a 8" - - |(% - A 4 ) 


, a 6 A9 
®9 = T " T 


a 9 " = - |(A 2 - A 5 ) 
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V =B 0 

V =B i 
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TABLE I.- CQKETXCISinS HJ THE EXPRESSIONS TOR 







! ° 6 " 

1 " 


-Q.04l67 


0.020855 0 
.020665 -.0000005 

.019642 -.0000002 | .000501 

000905 I .017425 -.0000652 I .000905 

.014465 -.0001772 

.011562 - .OOO2765 

.009015 --.0005109 
02877 | .000955 I .006905 -.0002959 


.005218 -.0002159 
.005894 -.0000997 


.50000 

.57558 

.64279 

.70711 

.76604 

.81915 

.06605 

.90651 

1.00000 



0.002085 

.002074 

.002028 

.001889 

.001701 

.001526 

.001552 

.001197 

.001065 

.000951 

.OQQ856 

.000779 


-0.01250 0 
-.01195 I .00096 
-.01032 
-.00751 
-.00540 
-.00395 
-.00298 
-.00227 
-.00179 
-.00143 
-.00117 
-.00100 
-.00087 
-.00076 
-.00069 
-.00054 



-0.06334 

-.08180 

-.07666 

-.06594 

-•05574 

-.04709 

-.04005 

-.05458 

-.02984 

-.02620 

-.02329 

-.02097 

-.01911 

-.01763 

-.01649 

-.01465 


.00009 

.00060 

.00226 

•004o6 

.00549 

.00643 

.00705 

.00739 

.00756 

.00761 

.00761 

.00757 

.00751 

.00745 

.00728 


0 
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-.00488 

-.01023 
-.01318 
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-.01396 
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-.01251 
-.01176 
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- .01047 .OOI74 
- .00989 .00139 

-.00945 .00146 
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TABLE I.- COEFITCIEHTS Kf THE EJtFKESBIORS TOR THE VELOCITY P 01 EHTIAL - Concluded 
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Po 
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|I|M 
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.07538 
.07844 

not. no 
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1 1 t 1 1 1 1 0 
'ui ro 'fS ^ i vn \j? ^ 5 ^ 
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(a) Q 


= a 0 ', bi', c 0 ', do', e 0 ', f 0 '. 


Figure 2.- Variation of the quantities Q as functions of Pq 
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Po 

(d) Q = ag', d^', e ^' , f 2 S * 3 ' > ^ 2 " ‘ 


Figure 2.- Continued. 
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Po 

(f) Q = a^', d ^' , e^' , a ^" , f^". 


Figure 2.- Continued 



2.- Continued 
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(z) Q — 1 , dy 1 , eg } dg , . 


Figure 2.- Continued 

















